Import the ESRI shapefile of German districts and the district attribute table. Join the two data frames, transform the CRS to “EPSG:3035” and check your changes.
You need to rename one of the id variables or adjust your join accordingly (“id = district_id”).
We want a first descriptive visual of the distribution of Covid-19 cases in Cologne and the surrounding districts. Calculate the number of Covid-19 cases in the last 7 days (cases_7days) by population (population) and multiply with 100,000 (in Germany usually called "7 Tages Inzidenzzahl).
Select Cologne (district_id == 5315), find the surrounding districts and plot the Covid Rate in Cologne and the surrounding districts.
You can use the dplyr function bind_rows to combine the two spatial objects “Cologne” and “Cologne Surroundings”.
Save your data set of Cologne and its surrounding districts as a ESRI Shapefile.